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Previous solutions for the ELF mode series are discussed briefly. Particular attention 
is paid to the height-gain functions. The excitation of the modes for vertical and hori- 
zontal dipole excitation is also considered. 



1. Introduction 

In a number of earlier papers [Wait, 1960a and b] 
the author discussed the theory of electromagnetic 
wave propagation in the extremely low frequency 
range (i.e., 1.0 to 3000 c/s). The model adopted is 
a smooth homogeneous earth surrounded by a con- 
centric reflecting layer which is to represent the 
ionosphere. The field was represented by a sum 
of modes for both vertical and horizontal electric 
dipole excitation. Essentially the same model has 
been used by Schumann [1954] but some of his results 
differ from those of Wait [1960a]. It is the purpose 
of this note to discuss certain extensions and limita- 
tions of the theory. 

Initially the earth is taken to be flat. This ide- 
alization simplifies the discussion. Furthermore, for 
ELF, influence of sphericity can be introduced later 
in the analysis. 

2. Formulation 

Choosing the usual cylindrical coordinate system 
(p,<t>,z) the earth's surface is at 2=0 and the iono- 
spheric reflecting layer is at z=h. Free space of 
electrical constant €o,Mo occupies the space 0<Cz<Ch. 
At z=0 and z=h, it is further assumed that the 
fields satisfy the Leontovich boundary conditions. 
Again this is not an essential restriction but it does 
simplify the discussion somewhat. Thus 
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in terms of the (complex) conductivities a g and a f of 
the earth and the ionosphere, respectively. The 
preceding forms for Z g and Z t are valid at ELF 



when the earth and the ionosphere are both equiv- 
alent to isotropic conductors [Wait, 1960b]. In the 
above, t}^12()tt ohms and A g and A t are dimension- 
less complex quantities. 

3. General Form of the Mode Series 

First the source is taken to be a vertical electric 
dipole (of current / and length ds) located at 2=z , 
where 0<2 <^- From a previous analysis the fields, 
in the space 0<^z<^h, are obtained from 
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where 
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is a reference field and W, S, and T are dimensionless 
factors which are expressible as a sum of modes. 
In the usual notation, the latter are given explicitly 

by 
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H" (2) and H[ 2) are the Hankel functions of the second 
kind of order zero and one. C n are the solutions of 
the mode equation 



and 



HziORiiOe-^^e-* 2 ™ (12) 

4. Height-Gain Functions 



While the above result was given in previous 
papers [Wait, 1960a and b] its significance was not fully 
discussed there. Furthermore, the heights of source 
and observer above the ground were taken to be zero. 
It is of interest to look into the nature of the solution 
under more general conditions. 

The function Wis considered first. It is the ratio 
of the actual vertical electric field E z to the reference 
field E of the same dipole source if it, and the ob- 
server, were on the surface of a flat perfectly con- 
ducting plane. If &p»l, corresponding to the 
"far-zone", it follows that 

RZ-M)* i (t4) ± S n Sre- i2 * s " p '%(z )fn(z) , (13) 

fl/\ n=0 

where X is the free-space wavelength. The vertical 
field E z is thus proportional to a sum of modes whose 
amplitudes are proportional to 
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At short distances, such that 

P « h ,W-,(l-i-^, (14) 

provided z=z —0. This is the expected behavior 
for a dipole if the reflecting layer (at z=h) was neg- 
lected. Numerical values of W were presented on 
previous occasions for the intermediate range where 
kp is unrestricted but z=z =0 [Wait and Carter, 
I960]. 

As can be seen from (4) , the height dependence of 
the individual modes is described entirely by the 
function f n (z). When s=0, we see that 
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For most applications to ELF propagation, |(7||» 
|A||, and thus/ n (0) can be replaced by unity. This 
approximation is contained implicitly in earlier work. 
For present purposes it appears to be preferable to 
define a "height-gain function" in the manner 
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Using (10) this can be rewritten in the form 
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it is seen that for the zero mode 
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For most ELF applications F (z) departs only slightly 
from unity. Also, it should be noted that if the 
ground was a perfect conductor the linear term in 
height is absent. Similar remarks apply to the 
higher modes. In fact, if A g =0, we have 

F n (z) = cos kC n z 
for all n. 

It is seen from (6) that the height-gain functions 
for the H^ component is identical to F n (z). In 
the case of the E p field it is desirable to define a height 
factor G n (z) as follows 
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For a given mode, this quantity is proportional to 
the ratio of the horizontal electric field at height z to 
the vertical electric field at height zero. From (11), 
it is seen that 
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For mode zero 
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Again, under most practical conditions, only the 
first term of this series need be retained. In fact, 
if z=0, we have exactly that 



G n (0) = A g 



(23) 



which is simply a consequence of the Leontovich 
boundary condition imposed at z=0. If the ground 
is perfectly conducting 



G n (z)=iC n sin kC n z. 
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5. Excitation of the Modes 

We now consider the excitation coefficient S„ 
defined by (7). Since 
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it readily follows that 
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If |A g |<c|A<|, which is the usual case, the second 
term in the denominator can be neglected. Then 
to the same approximation 
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where use has been made of ]i i = e 2ikCnh . Actually, 
(27) is exact if A^=0 and it can be regarded as an 
excellent approximation provided |Ag|<c|A*|. For 
practical applications at ELF, this may be further 
approximated to 



and 



5 =o since khC <^l 

d n ^l since 2khC n Q±2m. 

6. Further Extensions 



The extension to a vertical magnetic dipole source 
is very simple. In view of the symmetrical prop- 
erties of Maxwell's equations we can write down 
the corresponding field expressions. Thus 



H 2 =W*H 

H P =-S*H 
E,= T*H oV 



(28) 
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where IT*, S* and T* are dimensionless factors and 
H =i(llv\)Kds(e- ik >)lp (31) 

where K is the magnetic current in the vertical 
magnetic dipole of length ds. The explicit expres- 
sions for W*, S*, and T* are identical in form to 
those given by W, S, and T (i.e., (4) to (6)). Now, 
however, the modes are solutions of the equation 
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where m is an integer and 
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The superscript h is to 



AJ^— and A?^— . 
g A g A, 

denote the horizontally polarized nature of these 

reflection coefficients. 

In this case we note that 
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which is a very small quantity compared with unity. 
Consequently, the excitation of ELF waves by a 
vertical magnetic dipole on the ground plane is 
very poor. However, when the dipole is raised the 
situation is greatly improved. The height-gain 
function of the modes for this case is 
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Since A g is a small quantity, the second term of this 
series may be quite large even for small heights. 

The height -gain function F%[z) describes the 
behavior of both the H 2 and E<p components of the 
vertical magnetic dipole. The appropriate heigh t- 
gain function for the horizontal magnetic field 
component H p is defined by the function 
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in complete analogy to (21) for the electric dipole 
case. 

The excitation of the waveguide modes for a 
horizontal electric dipole was treated previously 
[Wait, 1960a]. For an electric dipole of moment 
Ids situated at z=z and oriented in the direction 
0—0, it was shown that 



E z = cos ct>E S 



(38) 



where E and S are given by (3b) and (5), respec- 
tively. Also 

rjoH 2 = -sin $E T* (39) 

where, explicitly, 
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For most applications the fields derived from S (i.e., 
the vertically polarized waves) are only significant. 
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The fields associated with T* (i.e., the horizontally 
polarized waves) are hardly excited at all when 
kz <g.l which is the usual case. 

The field of a horizontal magnetic dipole of 
moment Kds has the same structure as the fields of 
the horizontal electric dipole. Explicitly 



and 



H z =cos <I>H S* 

E z =sin <t>H T/r) . 
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The solution for the horizontal magnetic dipole has 
been discussed recently by Galejs [1961]. He 
obtains his results by an application of the reciprocity 
theorem as outlined by Wait [1960a]. 

The remarks in the preceding section concerning 
the height-gain functions for the vertical dipoles 
can be carried over to the horizontal dipoles without 
modification. 

The generalization of the foregoing results to 
include the effects of earth curvature is not difficult. 
Provided one stays away from the antipode of the 
source dipole it is only necessary to multiply the 
field expressions by 



\sm ej 



where 6 is the central angle subtended by the source 
and observer at the center of the earth. The justi- 
fication for this step was given previously. It might 



be added, however, that in addition to the restriction 
in the neighborhood of the antipode it is also neces- 
sary that 

Ci>h/a (43) 

and 
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It is not difficult to verify that these conditions 
are well justified at ELF. 
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